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dà 要 针对 穿孔 等 离子 弧 焊 接 工艺 , 建立 了 定点 焊 的 二 维 轴 对 称 非 稳 态 数 理 模型 , 描述 高 温 流动 的 电弧 与 工件 熔 池 的 夸 
合 输 运 过 程 , 通过 数值 模拟 获得 了 电场 、 人 磁场 \ 流 场 与 温度 场 的 演变 规律 . 模拟 焊 缝 熔 合 线 与 实验 结果 吻合 较 好 , 验证 了 数 
学 模型 的 合理 性 . 研究 结果 表明 : W 阴极 尖端 附近 的 电流 密度 与 温度 是 最 高 的 ; 等 离子 弧 在 工件 上 部 呈现 出 “ 钟 形 ”, 而 在 中 
心 轴线 处 因 进 入 小 孔 呈 细 长 的 “ 锥 形 ”; 高 速 Ar 气 等 离子 体 冲 击 到 小 孔 壁 面 时 速度 急剧 降低 , 小 孔 内 出 现 高 压 区 和 Ar 气 等 
离子 体 的 返 流 现象 ; 熔 池 内 流动 与 传 热 的 综合 作用 使 焊 颖 熔 合 线 呈 倒 “喇叭 形 ”. 步 对 焊接 的 操作 参数 及 焊 枪 结构 参数 
进行 了 影响 因素 的 正 交 试验 模拟 . 极 差分 析 表 明 , 焊 枪 的 结构 因素 比 焊 接 操作 参数 更 重要 , 即 对 电极 间距 、 电 极 内 缩 量 、 喷 
嘴 孔 径 需 予以 重点 考虑 , 以 获得 良好 的 焊 颖 形状 . 

关键 词 等 离子 弧 , 焊接 , 熔 池 , 耦合 , 正 交 试 验 
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ABSTRACT A 2D axial symmetrical mathematical model was developed for stationary keyholing plasma arc 
welding (PAW), to describe the transport process in coupled high-temperature flow arc and molten pool in the 
workpiece. The evolutions of electric, magnetic, velocity and temperature fields were simulated. The simulated fu- 
sion line of the weld bead is in quite good agreement with the experimental results, validating the mathematical 
model. It turns out that, both the current density and the temperature reach the maximum values near the tip of the 
tungsten cathode. The arc displays a typical bell-shape above the workpiece, but becomes slim cone-shape near the 
central axis as the arc enters the keyhole. The argon plasma slows down sharply when it strikes the inner wall of 


the keyhole, so high pressure appears in the keyhole and some argon plasma flows back. The combination of fluid 
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flow and heat transfer contributes to the reversed bugle shaped fusion line. The simulation of orthogonal test was 
further conducted to study the effects of operational and structural parameters of the weld torch. The range analysis 
shows that the structural parameters of weld torch are more influential than the operational parameters. That is, 


more attention should be paid to control the gap between two electrodes, the electrode shrinkage and the nozzle di- 
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ameter to guarantee the welding quality. 


KEY WORDS plasma arc, welding, molten pool, coupling, orthogonal test 


在 高 速 发 展 的 现代 加 工 制造 业 中 , 焊接 已 经 渗 
透 到 其 中 的 各 个 领域 . 随 着 科学 技术 的 进步 , 产品 
装备 对 轻 量化 、 节 能 化 、 高 参量 化 要 求 日 益 提高 , 这 
对 焊接 新 工艺 和 新 技术 提出 了 越 来 越 高 的 要 求 . 在 
当前 最 先进 的 3 种 高 能 束 流 " 焊 接 中 , 尽管 电子 束 焊 


追踪 了 定点 焊 小 孔 界 面 的 演变 过 程 , 模拟 实现 了 工 
件 的 穿孔 , 然而 该 模型 中 电弧 向 工件 传递 的 热量 以 
及 熔 池 中 气 - 液 界面 两 侧 气 相 与 液 相 的 流体 速度 都 
是 通过 既定 表达 式 给 定 的 , 没有 考虑 电弧 对 工件 炊 
池 的 热力 作用 . 张涛 等 "采用 烩 -孔隙 度 法 处 理 相 变 


与 激光 焊 在 精密 焊接 及 能 量 控 制 等 方面 具有 一 定 
的 优越 性 , 但 是 它们 设备 复杂 并 且 价 格 昂贵 , 同时 操 
作 要 求 较 高 , 不 适合 大 规模 焊接 生产 . 而 等 离子 弧 焊 
接 (plasma arc welding, PAW) 设 备 造价 及 运行 成 本 
较 低 , 焊接 的 操作 过 程 相对 简单 , 因而 被 人 们 极 大 
地 关注 . PAW 的 能 量 密度 可 达 10~10' W/cm’, 其 电 
弧 的 最 高 温度 可 达 52000 K”. 在 焊接 过 程 中 , PAW 
对 工件 的 变形 影响 小 , 工件 残余 应 力 小 , 焊接 成 型 
的 几率 比较 大 , 焊 颖 深 宽 比 大 尤其 是 小 孔 型 
PAW, 它 使 普通 PAW 方法 的 导热 焊 转 变 为 穿孔 焊 ， 
加 深 了 传 热 深 度 , 极 大 地 改善 了 焊 颖 附近 热量 分 布 
和 利用 , 能 够 形成 很 大 的 焊接 深度 和 很 小 的 热 影响 
K. 因此 , 小 孔 型 PAW 热量 集中 , 485338 7] 8, 焊接 速 
度 快 , 可 以 实现 单 面 焊 双 面 成 型 , 能 够 进行 10 mm 
以 下 壁 厚 材料 的 无 坡 口 对 焊 连 接 . 然而 , 小 孔 型 
PAW 也 有 许多 关键 问题 尚 待 解决 , 如 小 孔 对 焊接 工 
艺 条 件 的 变化 比较 敏感 ; 焊 缝 成 型 的 稳定 性 较 差 ; 
获得 良好 焊接 质量 的 合理 规范 参数 区 间 秆 等 后 . 
因此 , 近年 来 穿孔 PAW 相关 的 研究 工作 受到 了 广 
泛 重 视 . 
针对 电弧 部 分 , Tashiro 5:9! sr. f rs pi d 
(hollow cathode arc, HCA) 模 型 来 研究 不 同 工 况 下 
Ar 气 的 热源 特性 , 发 现 HCA 焊接 比 W 极 气体 保护 
电弧 焊 (gas tungsten arc welding, GTAW) 对 工件 的 加 
热 更 加 均匀 . Aithal 等 "建立 二 维 PAW 模型 , 模拟 等 
离子 体 在 焊 枪 内 的 层 流 流动 以 及 向 工件 表面 的 喷 
射 过 程 ,用 工件 表面 上 的 滞 止 压强 分 布 来 验证 模型 
的 可 靠 性 , 但 是 没有 考虑 热 辐射 . Tashiro 等 "研究 了 
当 弧 长 变化 时 焊接 的 电源 特性 对 弧 温度 与 弧 功 率 
的 影响 作用 . 

针对 PAW 燃 池 的 传 热 现象 , 人 们 开展 了 大 量 的 
研究 5, 但 往往 忽略 了 燃 池 的 流动 及 穿孔 . 为 研究 
炊 池 中 的 穿孔 行为 , 王 小 杰 等 利用 Level-Set 方 法 


潜 热 与 动量 损耗 问题 , 并 用 流体 体积 函数 (volume 
of fluid, VOF) 法 追踪 小 孔 界 面 . 霍 玉 双 等 "使 用 “ 双 
椭 球 体 + 峰值 热流 递增 的 圆柱 体 ” 的 组 合式 体积 热 
源 模型 来 确定 熔 池 的 形状 与 尺寸， 再 在 此 基础 上 通 
过 熔 池 内 部 的 受 力 分 析 建 立 小 孔 的 形状 方程 , 最 终 
确定 了 小 孔 的 形状 . IMR TESTE RI ODURSER E " 
组 合式 体积 热源 模型 以 及 焊接 熔 池 的 表面 变形 方 
T, 研究 小 孔 形状 随 受 控 脉 冲 电流 的 动态 变化 过 
T, 实现 了 小 孔 的 “打开 ”与 “闭合 ”Li 等 中 提出 复合 
体积 热源 模型 , 由 双 椭 球 热源 和 圆锥 体 热源 组 成 ， 
它 能 兼顾 热效应 和 小 孔 效 应 , 很 好 地 模拟 再 现 了 焊 
缝 熔 合 线 的 驼峰 . 但 以 上 研究 都 未 考虑 等 离子 弧 的 
实际 作用 , 是 用 简化 的 等 效 热 源 模型 来 描述 等 离子 
弧 提 供 的 能 量 和 穿孔 效应 . 

基于 等 离子 弧 与 燃 池 的 耦合 模拟 , Yamamoto 
等 I 研究 金属 蒸发 对 电弧 及 熔 池 的 影响 . Tashiro 
等 凹 建 立 了 考虑 熔 池 中 流动 的 弧 与 熔 池 的 耦合 模 
型 , 通过 预 设 圆 柱 形 小 孔 直 径 获 得 表面 张力 、 蔡 汽 
压力 、 切 应 力 随 小 孔 直 径 的 变化 规律 曲线 , 最 终 由 
三 力 平衡 确定 了 稳定 状态 下 的 小 也 直径 . Tanaka 
等 中 建立 了 二 维 轴 对 称 模型 , 研究 弧 与 熔 池 的 传 热 
流动 对 爆 颖 几何 的 影响 , 研究 表明 , 熔 池 的 流动 由 
电弧 拖 搜 力 、 电 磁力 、 浮 升力 、Marangoni 力 引起 ， 相 
kma, 电弧 拖 搜 力 与 Marangoni J Esx AEH, 二 
者 的 平衡 决定 了 熔 池 中 的 环流 方向 . 

综 上 所 述 , 尽管 对 PAW 的 模拟 研究 很 多 , 但 多 
是 对 工件 及 熔 池 单独 进行 的 , 无 法 真实 再 现 高 温 高 
速 等 离子 弧 对 熔 池 复杂 的 热力 作用 ; 在 弧 与 熔 池 耦 
合 的 研究 中 , 几乎 都 未 考虑 等 离子 弧 对 工件 的 穿 通 
孔 作用 , 也 未 考虑 熔 池 的 流动 , 影响 了 对 熔 池 中 流 
动 相 变 传 热机 理 的 准确 认识 . 因此 , 本 工作 拟 建立 
定点 PAW 中 等 离子 弧 与 熔 池 耦合 的 二 维 轴 对 称 数 
学 模型 , 分 析 电 弧 及 熔 池 中 的 电场 、 磁 场 、 流 场 、 温 
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度 场 的 演变 规律 及 相互 作用 ; 将 焊 颖 熔 合 线 模拟 结 
果 与 实验 数据 进行 对 比 , 以 验证 本 工作 模型 的 可 靠 


CD 处 喷 入 , 保护 气体 从 保护 气 喷嘴 的 入 口 GH ATE 
入 . 由 于 焊接 过 程 受 到 一 系列 物理 因素 的 复杂 作 


性 ; 并 进一步 对 焊接 的 结构 参数 、 操 作 参 数 进行 正 
交 试 验 设计 , 通过 模拟 和 极 差 分 析 确 定 各 影响 因素 
的 主 次 关系 , 以 期 为 实际 穿孔 PAW 工艺 优化 提供 
指导 . 
1 数理 模型 

图 1 是 穿孔 PAW 示意 图 . 等 离子 弧 焊 枪 的 W H 
极 和 工件 分 别 连 接 于 电源 的 两 极 , 通电 后 在 W 阴极 
和 工件 之 间 产 生 电 流 , 由 于 电流 的 Joule 热 作 用 使 
Ar 气 被 加 热电 离 形 成 等 离子 体 , 产生 了 电弧 . 电弧 
燃烧 于 W 阴极 和 工件 之 间 , 连续 送 入 的 工作 气体 穿 
过 电弧 空间 之 后 成 为 从 喷嘴 喷 出 的 高 温 等 离子 焰 
流 . 电弧 经 过 机 械 压缩 、 热 压缩 、 电 磁 收 缩 这 3 种 压 
缩 效 应 后 , 能 量 高 度 集中 在 直径 很 小 的 弧 柱 中 咏 . 弧 
柱 压缩 程度 较 强 , 等 离子 焰 流 喷 出 速度 较 大 , 因此 ， 
压缩 电弧 能 量 集中 , 穿 透 力 很 强 , 可 以 用 于 金属 材 
料 的 高 效 快速 焊接 . 
1.1 控制 方程 组 

穿孔 PAW 几何 模型 如 图 2 所 示 , AB 为 W 阴极 
Rim, 电极 内 缩 2 mm, 工作 气体 从 水 冷 喷嘴 的 入 口 


1 穿孔 等 离子 弧 焊接 (PAW) 示 意图 
Fig.1 Schematic of keyholing plasma arc welding (PAW) 


process 


图 2 穿孔 PAW 几何 模型 
Fig.2 Geometric model of keyholing PAW (unit: mm, r— 


radial coordinate, z—axial coordinate) 


用 , 为 突出 问题 本 质 , 建立 合理 的 模型 , 本 工作 引入 
以 下 假设 条 件 ™: (1) 等 离子 体 光 学 薄 , 且 处 于 局 
部 热力 学 平衡 状态 ; (2) 等 离子 体 的 热 导 率 、 密 度 、 
黏度 . 定 压 比 热 窜 和 电导 率 均 是 温度 的 函数 ; (3) 电 
弧 中 高 温 等 离子 体 与 炊 池 中 液态 金属 流动 均 为 层 
流 、 不 可 压缩 Newton 流体 ; (4) 预 设 一 个 贯穿 工件 
的 倒 圆 台 型 固定 小 孔 , 其 底部 半径 由 实验 数据 " 确 
定 ; 实验 发 现 焊 接 稳定 后 小 孔 底 部 的 尺寸 变化 较 
小 "2, 将 小 孔 底 部 半径 根据 实验 数据 取 为 定 值 
0.5 mm"; (5) 工件 中 除了 金属 的 热 导 率 以 及 液态 番 
度 与 温度 有 关外 , 其 它 热 物性 参数 均 为 常数 ; (6) X 
池 自 由 表面 为 平面 无 变形 , 忽略 其 蒸发 损失 . 

基于 以 上 假设 , 建立 了 如 下 的 等 离子 弧 与 炊 池 
输 运 过 程 的 二 维 轴 对 称 数 学 模型 . 


质量 守恒 方程 : 
9p 10 0 E 
P L L2 ps) (pv)=0 a) 
7 方向 动量 守恒 方程 : 
Xpv) 1909， ð - 
Pod eL gil) ps) 
1 0 Ov, a, Ov, OV h v, 
jg P rt à: V ðz T ar? da (2) 
z 方 向 动量 守恒 方程 : 
o(p.) 129 B. on 
-a "par Pee) tP) 
ĝ 9v.、 198 Qv, v. 
a” J rar” 0z 下 G) 
能 量 守恒 方程 
Xpc,T) 19 à - 
àt * zg PeT) F az PeT) = 
19 oT, 9 0T 
a a ài) t^ 4) 
电荷 守恒 方程 : 
18, op、0，op、 
erm Mrs irm Lig (5) 
Ohm 定律 
NE (6) 
or 
J.= ag e (T) 
0z 
做 矢量 势 方程 : 
198,94, A, E A, " 
7 ar ór )* az? B oJ, (8) 
19194 .04 
r Or 3 7 oz = Hod: (9) 
B-VxA (10) 


式 中 ,1 是 时 间 , s; p 是 密度 , kg/m’; rz 4) 200 2346 [a] 
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与 轴 向 坐标 , m; ww 和 分 别 为 径 向 与 轴 向 速度 , m/s; 
LÆ 7) 85 EE, kg/(m-s); % 和 中 分别 为 径 向 与 轴 疝 动 
量 守 恒 方程 源 项 , N/m c, 是 比热容 , (kg"K); TIE 
力学 温度 , K; k JE 5238, W/(m.K); 5. 是 能 量 守恒 方 
程 源 项 , Wim’; c 是 电导 率 , S/m; q 是 电势 ,V; J, MJ. 
分 别 为 径 向 与 轴 向 电流 密度 , A/m^; 4,. 和 4. 分 别 为 径 
问 与 轴 疝 磁 矢 量 势 ,Vs/m; mm 是 真空 磁 导 率 , H/m; B 
是 磁感应 强度 , T; A 是 磁 矢 量 势 , V -s/m. 

采用 答 - 多 和 孔 技术 名 处 理工 件 熔 化 过 程 的 相 变 . 
材料 总 烩 值 被 定义 为 显 烩 与 潜 热 的 总 和 , B: 
H=h+f, L, (11) 
式 中 , HÆRE AE, J/kg; h © TK, J/kg; L&K 
化 潜 热 , Jkg; f, 是 液体 体积 分 数 ，f 被 定义 为 : 


0 (T<T) 
Tef 
/= T-T. (T, «T«T) (12) 
1 (Tz T) 
XX, RÆ 304 不 锈 钢 的 固 相 线 温度 , K; T; 2 304 不 
锈 钢 的 液 相 线 温度 , K. 


两 相 区 ( 糊 状 区 ) 被 视 为 多 孔 介 质 , 考虑 固 / 液 相 


间 的 相互 作用 力 , 根据 Darcy 定律 "确定 : 
zak 
u= "m (13) 


AF, u 是 糊 状 区 流体 速度 , m/s; 及 是 渗透 率 ; 4 是 


液态 金属 动力 条 度 , kg/m.s); Vp 是 由 于 固 / 液 相间 
的 相互 作用 而 引起 的 压 降 , N/m’. 

应 用 Kozeny-Carman 表达 式 中 计算 该 区 域 的 渗 
IUE 


ENE À 
2 D. -fY ui 
式 中 , Do 是 与 多 孔 介 质 结构 有 关 的 系数 , 根据 文献 
[8.25], 模型 中 可 取 为 常数 1.2x10". 
考虑 式 (12)~(14), 在 熔 池 动量 守恒 方程 中 加 入 
如 下 源 项 , 以 表征 固 / 液 相 界面 上 的 相间 作用 力 : 


= 
S= K” (15) 


式 中 , S$ 是 固 / 液 相 界面 上 的 相间 作用 力 , N/m’. 

动量 与 能 量 守 恒 方 程 中 的 源 项 表达 式 见 表 1. 
电弧 动量 守恒 方程 源 项 中 主要 有 电磁 力 , 能 量 守 恒 
方程 的 源 项 有 Joule 2. E Y Aja 4; 868 8] 38 R. 
磁 矢 量 势 微分 控制 方程 及 电荷 守恒 方程 的 求解 正 
确 , 是 电弧 中 动量 守恒 方程 中 电磁 力 与 能 量 守 恒 方 
旦 中 Joule 热 正确 加 载 的 保证 与 前 提 . 磁场 的 求解 采 
j 磁 矢量 势 微分 控制 方程 , 避免 了 按照 Ampere 环 
路 定理 对 电流 密度 进行 复杂 的 积分 计算 磁场 . V 
动量 守恒 方程 中 考虑 切 应 力 、 重 力 和 Darcy 阻力 , 忽 
KEREJ”. 

H A R R ih S R27], 如 式 (16) 


>H 


Uml 


R1 动量 与 能 量 守恒 方程 的 源 项 


Table 1 Source items in momentum and energy equations 


Zone Source item Expression 
ya 
Ar S, posue em 
e ' ðr J.B, 
S E pgtT k 
z Oz 8 tk 
2 72 
S J, tJ; Ske oT | 0T 
° Cr 2e Ua Joz) Á 
op pH 
Molt ] S Earl 
olten poo. » y K^" 
op 从 
S UP pog 
Z az PE- K” 
S H &(pf,L.) 


-lð -3 
a - 13 ipo, fL)- So fiL) 


Note: p—pressure, Pa; r, z—radial and axial coordinate, respectively, m; J,, /一 radial and axial compo- 


nent of the current density respectively, A/m^; Bs—the azimuthal component of magnetic field, T; 


p—the density, kg/m^; e—the acceleration due to gravity, m/s; o—the electrical conductivity, S/m; 


ks—Boltzmann constant, J/K; e—electronic charge, C; 7—thermodynamic temperature, K; R—radi- 


ation loss per unit volume, W/m^; uj—dynamic viscosity of liquid metal, kg/(m*s); K—permeabili- 


ty; v, v.—radial and axial velocity, respectively, m/s; f.—liquid fraction; L,—latent heat of fusion, 


J/kg; t—time, s 


一 个 一 
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dame - Avar 

R-4.17162 x 10?T^ - 2.33445 x 107T^ + l 2m ór (19) 
0.004397 - 17.00853 (16) - LL 


1.) 定 解 条 件 
本 工作 计算 域 如 图 2 所 示 , 其 中 4BCDEF- 
GHIMKLA 为 等 离子 体 区 域 , JKMI 为 工件 区 域 . 模 
型 计算 条 件 如 表 2 所 示 . 其 中 , n 表示 法 向 方向 ; LA 
是 对 称 轴 ; AB 是 焊 枪 的 电流 入 口 ; BC 是 焊 枪 电极 壁 
H; DE, EF 和 FG 是 水 冷 Cu 壁 面 ; KL FUR KE H 
O; HT 是 上 部 等 离子 体 出 口 ; CD 和 GH 分 别 是 工作 
气体 和 保护 气体 的 速度 入 口 , 且 假 设 为 等 速 分 布 ; 
和 JK 分 别 是 工件 侧 表 面 和 下 表面 , 其 热 边界 条 件 
如 式 (17) 所 示 , 考虑 了 对 流 换 热 与 辐射 换 热 : 
q7- (T, - T.)* ea(T, - T.) (17) 


式 中 , q 是 热流 密度 , W/m2 e3& Jb iE EE TR 8 I 3E; a 
是 Stefan-Boltzmann 常数 , W/(m^-K^; T LFR h 
温度 , K; T, 是 环境 温度 , K; h 是 表面 对 流 换 热 系 数 ， 
W/(m"-K). T.73300 K, hh, 的 计算 参考 文献 [28]: 
h,22.169* 10 ^7^* (18) 


MI 和 MK 分 别 为 工件 上 表面 和 小 孔 侧 壁面 . 由 
于 焊接 熔 池 表面 温度 分 布 不 均匀 , 表面 张力 随 温度 
梯度 变化 产生 Marangoni 力 站 ( 式 (19)), 它 对 熔 池 的 
流动 驱动 起 重要 作用 


表 2 初始 条 件 和 边界 条 件 


式 中 , s 和 二分 别 为 径 向 与 轴 向 切 应 力 , N/m’; y 是 表 
面 张力 , N/m, 它 的 数值 由 实验 数据 中 给 定 . 
2 数值 求解 

本 工作 采用 计算 流体 力学 (CFD) 软 件 FLUENT 
对 数学 模型 进行 求解 . 采用 非 结 构 化 网 格 对 计算 区 
域 进行 离散 化 , 靠近 对 称 面 处 选用 小 尺寸 密集 网 
格 , ren DX el ide FEL STE PL e. 采用 PISO 算法 
求解 离散 的 微分 控制 方程 组 , 其 中 动量 守恒 方程 和 
能 量 守 恒 方 程 的 对 流 项 采用 二 阶 迎 风 离 散 格 式 , 压 
力 插值 采用 Standard 格式. 数值 计算 收敛 条 件 是 连 
续 性 方程 的 绝对 残 差 和 小 于 107, 而 能 量 方程 绝对 
残 差 和 小 于 10 其 余 方 程 的 绝对 残 差 和 均 为 小 于 
10 同时 , 为 避免 假 收 敛 现象 , 还 补充 MT 和 MK 界 
面 上 两 侧 总 热流 量 绝 对 值 的 相对 差 值 小 于 10” 作 为 
能 量 方程 的 收敛 判 据 . 时 间 步 长 设置 为 0.01 s. 
民 据 现 有 研究 ,等 离子 弧 物性 参数 受 温度 
影响 很 大 , 本 工作 采用 文献 34] 给 出 的 Ar 气 等 离子 
体 密度 .比热容 .黏度 、 热 导 率 及 电导 率 随 温度 的 变 
化 值 , 如 图 3 所 示 . 表 3 给 出 了 304 不 锈 钢 平 板 的 物 
性 参数 与 焊接 工艺 参数 的 数值 . 


Table 2 Initial and boundary conditions 


Condition y, y. T 9 A 
Initial 0 0 300K 0 0 
Boundary 
AB 0 0 3000 K -24 V ðA/ðn=0 
BC 0 0 3000 K ðg/ðn=0 ðA/ðn=0 
CD 0 0.758 m/s 300 K ðp/ðn=0 ðA/ðn=0 
DE 0 0 1000 K ðp/ðn=0 0A/àn 2 0 
EF 0 0 1000K ðp/ðn=0 0A/àn 2 0 
FG 0 0 1000K ðp/ð3n=0 ðA/ð3n=0 
GH 0 0.764 m/s 300 K ðg/ðn=0 0A/àn 2 0 
HI àv, [0n 2 0 0 300K ðp/ðn=0 ðA/ðn=0 
IJ 0 0 Eq.(17) 0 0 
JK 0 0 Eq.(17) 0 0 
KL 0 ðv./ðn=0 300 K op/on=0 0 
LA óv,/ón - 0 àv. /ón 2 0 9T /àn 20 ðg/ðn=0 9A/àn - 0 


Note: LA—axis of symmetry; AB—current entry of welding torch; BC—wall of the electrode; DE, EF, FG— 


wall of water-cooled Cu; KL—bottom outlet of keyhole; H/—upper exit of plasma; CD, GH—entry 


of working gas and protective gas; JJ, JK—side surface and lower surface of workpiece; o—electric 


potential, V; A—magnetic vector potential, V *s/m; n—normal direction 
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Temperature / 10? K 
3 Ar 气 等 离子 体 的 热 物性 随 温度 变化 曲线 
Fig.3 Temperature-dependent thermophyical properties of Ar plasma" 
(a) density and dynamic viscosity 
(b) specific heat 
(c) thermal conductivity and electrical conductivity 
表 3 304 不 锈 钢 平板 的 物性 参数 与 焊接 工艺 参数 w 
Table 3 Thermophyical properties of stainless steel 304 and welding operational parameters"! 
Nomenclature Symbol Value 
Solidus temperature T. 1663 K 
Liquidus temperature Ti 1723 K 
Latent heat of fusion L, 2.6x10? J.kg” 
Specific heat Cimer 630 J-kg 天” 
Density pom 7200 kg-m^? 
Elementary charge e 1.602x10^ C 
Boltzmann constant ks 1.381x10? J-K“ 
Magnetic permeability Lo 4nx10" H.m 
Radiation emissivity E 0.4 
Welding current i 169 A 
Arc voltage u 24 V 
Working gas flow rate qw 3 L-min” 
Shielding gas flow rate qs 30 L-min ' 
Gap between two electrodes E^ 5mm 
Electrode shrinkage k 2 mm 
Nozzle diameter d, 2.8 mm 
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OČE). 这 些 作用 使 得 电弧 能 够 达到 极 高 温度 . W 
阴极 尖端 处 电流 密度 最 高 , 达到 9.9x10* A/m’. 沿 径 


实际 工艺 和 模拟 结果 均 表明 电弧 温度 在 整个 向 偏离 弧 柱 中 心 轴 的 区 域 电流 密度 大 幅度 减 小 , 这 


过 程 中 变化 很 小 , 即 在 穿孔 等 离子 弧 定 点 焊接 过 程 。 是 由 于 喷嘴 对 高 温 缴 柱 实施 的 机 械 压缩 和 热 压 缩 


中 电 最 容易 维持 热 稳定 状态 . 图 4-7 给 出 了 2.0 s 时 作用 使 得 该 区 域 温度 较 低 , 因而 该 区 域 Ar 气 等 离子 
刻 弧 柱 及 工件 中 各 物理 量 的 分 布 状态 , 包括 电场 、 。 体 电导 率 很 低 电流 几乎 无 法 从 该 区 域 通过 . 也 下 


电磁 力 场 ,温度 场 ` 流 场 与 压力 场 . 


分 布 . 焊接 时 所 加 载 的 电压 为 24 V, 工件 作为 阳极 ， 
由 于 接地 而 成 为 零 等 势 体 . 关注 焊 枪 喷嘴 中 W 阴极 


图 4 显示 了 


是 此 作用 使 得 电弧 的 电流 积聚 在 喷嘴 内 部 , 从 而 拥 
高 能 量 密度 来 实现 穿孔 作业 . 

图 5 显示 了 弧 柱 中 电磁 力 的 分 布 . 从 电磁 力 的 方 
向 可 以 看 出 , 电磁 力 对 电弧 起 收缩 与 加 速 的 作用 . 


弧 柱 中 的 电势 (无) 与 电流 密度 ( 右 ) 


附近 区 域 : 由 图 6a 和 b 可 知 , Ar 气 等 离子 体温 度 很 W 阴极 尖端 附近 电磁 力 最 大 , 达到 2.6x10 N/m. 从 
高 , 则 此 处 Ar 气 等 离子 体 的 电导 率 也 很 大 , 参考 图 。 所 量 图 中 看 出 , 电磁 力 对 电弧 主要 起 收缩 作用 , 使 
3c; 同时 从 图 4 中 看 出 , 该 区 域 电势 梯度 很 大 , 因而 ” 弧 柱 直径 变 细 , 提高 其 能 量 密度 . 

形成 强 电 流 密度 , 进而 产生 大 量 的 Joule 热 ( 即 图 6a 是 电弧 及 工件 中 速度 ( 左 ) 与 温度 ( 右 ) 分 
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4 在 2.0s 时 刻 弧 柱 及 工件 中 电势 与 电流 密度 分 布 


Fig.4 Distributions of electrical potential and current density in arc column and workpiece at 2.0 s 


16 
Maximum E d 14 
2. 6*10' Nim 7 pp 


图 $ 在 2.0s 时 刻 弧 柱 中 电磁 力 分 布 


Fig.5 Distribution of electromagnetic force in arc column at 2.0 s 
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6 在 2.0s 时 刻 弧 柱 及 工件 中 速度 与 温度 分 布 


Fig.6 Distributions of velocity and temperature in arc column and workpiece at 2.0 s (T;.—temperature of arc, Ts—tem- 


perature of workpiece) 

(a) whole calculation domain 

(b) local region close to cathode 

(c) distribution of velocity upon the top surface of workpiece 


(d) distribution of velocity in the molten pool 


fg. 对 于 弧 柱 , W 阴极 尖端 附近 的 温度 最 高 , 达到 ”上 典 的 “ 钟 形 ”, 而 在 中 心 轴线 处 因 进 入 小 孔 呈 细 长 的 
35000 K, 如 图 6b 所 示 . 喷嘴 的 压缩 作用 以 及 保护 气 “ 锥 形 ”. 喷嘴 出 口 处 的 弧 柱 中 心 轴 线 上 轴 向 速度 最 
的 作用 , 使 弧 柱 等 温 线 在 工件 上 部 的 外 围 旦 现 出 经 。 大 , 达到 850 m/s, 与 文献 [35] 中 电弧 速度 相符 . 高 温 
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7 在 2.0s 时 刻 弧 柱 9 


P 流 线 与 压强 分 布 


Fig.7 Distributions of velocity streamline and pressure in arc column at 2.0 s (p—pressure) 


与 高 速 的 结合 使 弧 柱 中 温度 分 布 更 加 均匀 , 强化 等 
离子 弧 射 流 的 焊接 效率 与 质量 . 模拟 结果 表明 , 18 
弧 柱 中 心 轴线 方向 在 离开 W 阴极 尖端 3~12 mm 范 
EA, 电弧 温度 梯度 很 小 , 也 就 是 说 从 喷嘴 出 口 到 
小 孔 内 部 Ar 气 等 离子 体 的 温度 几乎 保持 一 样 高 , 有 
利于 热量 沿 金属 厚度 方向 传输 , 提高 其 炊 透 速率 . 
少 部 分 Ar 气 等 离子 体 从 小 孔 底 部 流出 , 温度 和 速度 
都 很 高 , 速度 几乎 可 达 300 m/s; 大 部 分 Ar 气 等 离子 
体 从 工件 上 表面 流出 , 如 图 6c 所 示 , 其 温度 和 速度 
都 相对 较 低 . 高 温 等 离子 弧 通 过 导热 、 对 流 与 辐射 
方式 将 热量 传递 至 工件 . 当 工 件 温度 高 于 金属 液 相 
线 温度 ( 即 熔 化 温度 , 1723 K), 工件 熔化 形成 熔 池 . 
燃 池 在 Marangoni 力 及 切 应 力 的 作用 下 产生 环流 流 
动 , 如 图 6d 所 示 , 最 大 流动 速度 达到 0.15 m/s. Ke 
受 对 流传 热 的 影响 , 上 部 熔 宽 增 大 , 熔化 轮廓 线 ( 熔 


冷 喷嘴 及 保护 气 对 其 压缩 作用 尚未 表现 出 来 . 0.4 s 
时 , 高 温 电 弧 在 工件 上 表面 被 压低 ,“ 钟 形 ” 的 电弧 
轮廓 形成 . 0.6-0.8 s EF, 弧 柱 中 外 围 电弧 进一步 被 压 
低 , 并 出 现 尖 端 , 尖端 与 弧 柱 中 轴线 的 距离 为 10 mm 
左右 .水冷 喷嘴 下 方 的 低温 区 增 大 , 温度 越 低 , Ar ^ 
等 离子 体 热 导 率 越 小 , 低热 导 率 阻碍 了 热量 的 外 
散 , 保证 了 弧 柱 中 心 附近 的 高 温 . 1.0 s 以 后 , 弧 柱 温 
度 场 基 本 稳定 . 弧 柱 中 的 速度 场 一 直 在 变化 . 高 速 
Ar 气 等 离子 体 喷射 到 小 孔 内 部 斜 壁面 上 时 , 速度 
迅速 减 小 , 出 现 Ar 气 等 离子 体 返 流 . 由 于 速度 急剧 
降低 , 在 小 孔 内 部 会 出 现 部 分 高 压 区 , 阻碍 Ar 气 等 
离子 体 的 向 下 冲击 , 700 m/s 以 上 速度 区 域 沿 轴 向 逐 
渐 向 喷嘴 回 缩 . 高 速 等 离子 弧 的 直径 收缩 得 非常 
细 , 1.0 s 以 后 , 300 m/s 以 上 的 高 速 等 离子 弧 都 集中 
在 直径 不 到 2 mm 的 柱状 区 域内 , 使 弧 柱 的 能 量 密 


合 线 ) 呈 倒 “ 喇 叭 形 ”. 工件 中 未 熔化 区 域 传 热 只 有 导 
热 作用 . 

图 7 是 弧 柱 中 流 线 ( 左 ) 与 压强 ( 右 ) 分 布 . 高 速 Ar 
气 等 离子 体 冲 击 到 小 孔 内 壁面 时 速度 急剧 降低 , 出 
现 部 分 高 压 区 和 Ar 气 等 离子 体 的 返 流 现象 . 喷嘴 出 
口 下 方 处 极 高 速 的 Ar 气 等 离子 体会 形成 对 保护 气 的 
卷 吸 作用 , 从 而 保护 气 可 以 压制 热 等 离子 体 而 保护 
工件 , 并 防止 高 温 熔 化 工件 吸收 空气 中 有 害 气 体 . 
3.2 弧 柱 与 工件 内 流 场 \ 温 度 场 的 演变 过 程 

图 8 显示 了 电弧 与 工件 中 速度 场 ( 左 ) 与 温度 场 
( 右 ) 随 时 间 的 演变 过 程 , 不 同时 刻 等 温 线 的 条 数 及 
其 上 对 应 的 温度 值 都 是 相同 的 , 这 样 有 利于 清晰 地 
比较 分 析 温 度 场 随 时 间 的 变化 情况 . 

在 电弧 区 , 初始 阶段 0.2s 时 电弧 自由 凸 起 , 水 


度 更 加 集中 . 

在 工件 区 , 0.6s 时 工件 开始 熔化 , BURG 2 
( 即 为 1723 K 的 液 相 线 ). 1.0s 时 , 熔 池 变 大 , 由 于 
熔 池 流动 较 弱 , 流动 对 燃 池 内 传 热 影响 较 小 . 1.4 s 
时 刻 , 熔 池 中 出 现 明显 的 环流 流动 , 这 是 由 于 Ma- 
rangoni 力 的 作用 引起 的 . 随后 熔 池 逐渐 变 大 , 同时 
热 影响 区 也 加 深 ; 熔 池 内 流动 速度 也 逐渐 增 大 , 最 
大 流速 达到 0.15 m/s. 1.8 s 之 后 , 整个 工件 被 完全 
AE, 熔 合 线 到 达 小 孔 的 下 表面 , 熔 池 内 流动 加 强 
并 促进 了 熔 池 中 的 传 热 , 出 现 明显 的 倒 “ 喇 叭 形 ” 
RER. 
3.3 实验 验证 

为 了 验证 模型 的 可 靠 性 , 本 工作 对 模拟 所 选用 
的 厚度 为 6 mm 的 304 不锈钢 平板 进行 了 定点 穿孔 
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8 弧 柱 与 工件 中 的 速度 场 、 温 度 场 动态 演变 过 程 


Fig.8 Dynamic evolution of velocity and temperature fields in the arc column and workpiece (t—time) 


(02s  (bO4s (c)O6s (d)O8s (e)yLOs 


PAW 实验 . 实验 中 焊接 工艺 参数 与 计算 模型 参数 
一 致 , 见 表 3™. 

焊接 结束 后 切 开 焊 件 , 通过 焊 缝 横 截 面 的 宏观 
金 相 照片 , 获得 了 焊 颖 的 实验 几何 形状 与 尺寸 , 与 
计算 结果 对 比如 图 9 所 示 . 可 见 , 2 者 形状 趋势 比较 
一 臻 . 上 部 炊 宽 相对 误差 为 1.22%, 底部 熔 宽 相对 误 
差 为 4.72%. 可 见 , 模拟 计算 与 实验 结果 吻合 良好， 
验证 了 本 工作 模型 的 正确 性 . 


(12s (gi4s (h)l.6s (1.8s 0)20s 


4 模拟 正 交 试验 分 析 

为 了 对 焊 枪 设计 以 及 PAW 焊接 工艺 提供 优化 
指导 , 本 工作 对 其 中 的 结构 参数 与 工艺 操作 参数 进 
行 模拟 正 交 试验 分 析 . 结构 参数 包括 焊 枪 的 电极 间 
距 、 电 极 内 缩 量 、 喷 嘴 孔 径 以 及 小 孔 的 上 表面 半径 ， 
操作 参数 有 焊接 电流 、 工 作 气 体 流量 、 保 护 气流 量 . 
总 共 7 个 因素 , 每 个 因素 3 水 平 , 如 表 4 所 示 . 对 于 3 
水 平 7 因素 的 正 交 模拟 分 析 , 需要 设计 18 组 试验 ， 
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9 焊接 熔 池 横断 面 实验 与 计算 结果 对 比 图 


Fig.9 Comparison of weld bead cross-section between experimental image (left) and calculated result (right) (The work- 


piece is stainless steel 304 and its thickness is 6 mm; the welding current is 169 A; the arc voltage is 24 V; the work- 


ing gas flow rate is 3 L/min; the shielding gas flow rate is 30 L/min; the gap between two electrodes is 5 mm; the 


electrode shrinkage is 2 mm; the nozzle diameter is 2.8 mm) 


表 4 正 交 试 验 模拟 因素 水 平 表 


Table 4 Table of factors and levels for orthogonal test simulation 


Level i/A q«/(L-min?) q./(L-min" g/mm k / mm d, / mm r, / mm 
1 169 3 20 5 2 2.8 2 
2 159 1 10 3 0 1.8 1 
3 179 5 30 T 4 3.8 3 


Note: i—current; gv 一 working gas flow rate; q.—shielding gas flow rate; g—gap between two electrodes; 


/一 electrode shrinkage; d,—nozzle diameter; r,—keyhole top radius 


记 为 Lis(3”. 由 于 在 实际 的 焊接 实验 9 
宽度 的 观测 比较 困难 , 这 使 得 其 实际 大 小 无 法 被 准 


H. ILEX É 


确 确定 , 因此 通过 模拟 来 预测 小 孔 上 表面 尺寸 的 可 


能 大 小 . 


为 了 研究 PAW 焊接 各 因素 的 主 次 关系 , 对 试验 


中 的 7 个 因素 进行 极 差 分 析 , 如 图 10 所 示 . 以 爆 颖 


的 上 部 与 底部 炊 宽 作为 考核 指标 , 因为 它们 最 能 


映 焊 颖 的 几何 形状 及 焊接 质量 . 从 图 中 可 以 看 出 ， 
TE PAW 焊接 过 程 中 , 电极 间距 、 电 极 内 缩 量 ` 喷 路 孔 
径 对 应 的 上 、 底 部 将 宽 极 差 数 是 较 大 的 ; 即 焊 枪 的 


结构 参数 对 焊 缝 熔 宽 影响 很 大 ， 


它们 应 该 是 PAW 焊 


接 工艺 的 重点 优化 内 容 . 尤其 是 , 对 应 图 中 最 大 极 
焊接 过 程 中 必须 慎重 确 
中 电流 密度 分 布 及 对 弧 


差 数 的 电极 内 缩 量 , 在 PAW 
E, 因为 电极 内 缩 量 对 弧 柱 


柱 压 缩影 响 其 大 , 进而 影响 


高 温 电弧 向 工件 的 冲击 


与 传 热 , 最 终 显 著 影响 焊 颖 的 熔 宽 . 相 比 于 结构 参 


数 , 焊接 电流 、 工 作 气 体 流 革 


ER CUL 


等 工艺 操 


作 参 数 对 应 的 上 、 底 部 熔 宽 的 极 差 数 较 小 , 即 焊 颖 


5 
—4— Top width range 
—e— Bottom width ra 


Range / mm 
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10 各 影响 因素 下 的 焊 颖 上 、 底 部 熔 宽 的 极 差 图 


Fig.10 Range analysis of top and bottom weld bead widths 


under different factors 


熔 宽 受 这 些 工 艺 操作 参数 的 
5 结论 


Ei^ 
HZ 


响 小 于 焊 枪 的 结构 


(1) 等 离子 弧 的 外 围 在 工件 上 部 呈现 出 经 典 “ 钟 
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E". 而 在 中 心 轴 线 处 因 进入 小 孔 呈 细 长 的 “ 锥 形 ” 
W 阴极 尖端 附近 Ar 气 等 离子 体能 量 密度 高 度 集 中 ， 
该 处 温度 最 高 可 达 35000 K; 等 离子 弧 的 高 速 区 集 
P 在 弧 柱 的 中 心 轴线 上 ; 当 高 速 Ar 气 等 离子 体 冲击 
到 小 孔 壁 面 时 , 其 速度 急剧 降低 , 小 孔 内 出 现 高 压 
区 和 Ar 气 等 离子 体 的 返 流 现象 . 

(2) 熔 池 在 Marangoni 力 及 切 应 力 等 作用 下 产 
生 环 流 流动 , 焊 颖 炊 合 线 呈 倒 “ 喇 叭 形状. 

(3) 相对 于 工艺 操作 参数 , 焊 枪 的 结构 参数 对 焊 
颖 熔 宽 (反映 焊 颖 几何 形状 及 焊接 质量 ) 的 影响 更 显 
著 , 需要 予以 重点 考虑 电极 间距 、 电 极 内 缩 量 和 喷 
嘴 孔 径 . 
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